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Drag resistance of 2D electronic microemulsions 

Boris Spivak^y and Steven A. Kivelson 2 'Lj] 
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Motivated by recent experiments of Pillarisetty et al, Phys. Rev. Lett. 90, 226801 (2003), we 
present a theory of drag in electronic double layers at low electron concentration. We show that 
the drag effect in such systems is anomolously large, it has unusual temperature and magnetic field 
dependences accociated with the Pomeranchuk effect, and does not vanish at zero temperature. 



Over the last decade a large number of dramatic phe- 
nomena have been discovered in studies of the two dimen- 
sional electron gas (2DEG) in MOSFET JJQ and semi- 
conductor heterojunction devices [3, OLE!- Collectively, 
these phenomena are referred to under the euphonious ti- 
tle of the "2D apparent metal insulator transition." They 
are associated with the behavior of the 2DEG in very 
clean (high mobility) devices with strong interactions, 
i.e. when the ratio of the typical interaction strength to 
the Fermi energy, r s = l/y/nog, is large, r s ~ 10 — 40. 
(Here n is the electron density per unit area and ap is the 
effective Bohr radius.) These observations make it clear 
that the physics of this simple, but conceptually vital 
system is much richer than was previously appreciated. 
Elsewhere, we have proposed 0, Q that these phenom- 
ena are associated with a set of "electronic microemul- 
sion phases" which occur as a consequence of Coulomb 
frustrated phase separation 9j in the ideal (zero disorder) 
2DEG at densities intermediate between the high densi- 
ties, where the system is a Fermi liquid, and low densities, 
where it is a Wigner crystal. 

In the present paper we concentrate on the drag resis- 
tance of 2D large r s metals, with the goal of explaining 
the unexpectedly large value of the drag resistance and 
the unusual temperature and magnetic field dependances 
that have recently been reported in bilayer p-type GaAs 
heterojunction devices [iaul- In a bilayer device, in the 
presence of a dissipative current I a in the first (active) 
layer, the inter-layer electron-electron interactions induce 
a momentum transfer between layers, which in turn pro- 
duces a voltage, Vp, across the second (passive) layer. 
In the linear regime this effect can be characterized by a 
drag resistance (per square) 



Rd = Vp/Ia 



(1) 



Measurements of Rd are a potentially useful probe of 
interesting correlations of an electron fluid. It is exceed- 
ingly small at low temperatures in a Fermi liquid while, 
as we shall see, it is parametrically larger in a "bubble 
fluid" phase. It is therefore an ideal diagnostic for the 
presence of such phases. 



The existing theories of drag [1{J, [HI are mostly based 
on the Fermi liquid theory of the electron liquid in indi- 
vidual layers. A common feature is that, at small tem- 
peratures, Rd is small and vanishes as T — > or as the 
spacing between layers, d — > oo: 



Rd ~ {h/e 2 )(k F d)- a *(T/E F ) a T. 



(2) 
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For a Fermi liquid, ap — 2 up to logarithmic corrections 
and ad = 2 — 4 depending on the value of kp£, where £ is 
the elastic mean-free path and kp is the Fermi momen- 
tum. 

Experimental Background: The drag experiments 
that motivated this study are performed on high mobil- 
ity samples whose conductances are significantly larger 
than e 2 /%. Therefore, the putative electron localization 
length one would infer on the basis of localization theory 
is exponentially long, and hence the associated 2D local- 
ization phenomena can be ignored for present purposes. 

In double layer devices with relatively small r s , 
measurements |12| of the drag are in qualitative agree- 
ment with Fermi liquid theory |lCtlll| . However, at large 
r s , experiments on p-GaAs double layers 0,13 differ sig- 
nificantly from the predictions of this theory : 

1. The drag resistance in these samples found to be 
2-3 orders of magnitude larger than expected on the basis 
of Fermi liquid theory. 

2. In the presence of a magnetic field, H\\, parallel to 
the film, the temperature dependence of Rd (T) is signif- 
icantly suppressed. (Here, H» primarily couples to the 
spins of the electrons, rather than to their orbital mo- 
tion. ) Whereas in a Fermi liquid, Rd{T) is a quadratic 
function of T, in large r s devices, at Hu = 0, the value 
of ap in Eq. [21 appears to be noticeably larger than 2; 
for example, or = 2.7 in [a, 111- The value of ap(H\\) 
then decreases with H\\ and saturates for H\\ > H* at a 
value which is significantly smaller than 2; for example, 
in Refs. @ and [||, a T (-H|| > H*) ~ 1.2. H* is gener- 
ally believed to be the field required to fully polarize the 
electron spins. 

3. As a function of increasing H\\, Rd(H») increases 
by a factor of 10-20 and saturates when H\\ > H* (See 
Fig. 1 in [5j] ) . In the framework of the Fermi liquid theory, 
one would expect a significant decrease of Rd{H\\) in the 
spin polarized state because Ep is increased by a factor 
of 2 as H\\ increases from to H\\ > H*, thus decreasing 
the electron-electron scattering rate. 



4. The T and especially the Hu dependences of 
Rd{T,H\\) and the resistances of the individual layers 
R(H\\ , T) look qualitatively similar to one another. (See 
Figs.l a,b in ) 

In fact, R(Hu,T), itself, appears highly anomalous 
from the viewpoint of Fermi liquid theory. Thus, al- 
though our primary focus will be on Rp>, we will also 
summarize the most salient anomalies in R(Hu , T). Since 
similar phenomena have been observed in studies of the 
2DEG in MOSFETs and of the 2D hole gas (2DHG) in 
semiconductor heterojunction devices, we report obser- 
vations from both systems: 

5. An apparent continuous metal insulator transition 
as a function of n has been observed in Si MOSFET's 

l| and to a lesser extend in GaAs heterojunctions 0, 
J. (No such transition is predicted on the basis of 
Fermi liquid theory.) For the most part, in this paper we 
will focus on the slightly higher density samples on the 
metallic side of the apparent transition. 

6. The resistance of these samples increases sig- 
nificantly with increasing T: by as much as a factor 
of 6 in Si MOSFET's Q and as much as 3 in GaAs 
heterojunctions 0, y| ■ 

7. The resistance of metallic samples at low T is sig- 
nificantly increased by non-zero Hu , and becomes nearly 
H\\ independent at high magnetic fields, Hu > H* pj. 
Moreover, non-zero Hu strongly suppresses the T depen- 
dence of i?, so much so that R is nearly temperature 
independent for H\\ > H* . (To date, this latter effect 
has been documented only in Si MOSFET's Q, where 
dR/dT at low temperature is decreased for H\\ > H* by 
as much as a factor of 100 relative to Hu = 0.) 

Theoretical Background: The phase diagram of the 
bilayer 2DEG depends on many parameters. For sim- 
plicity, in this article we consider only the case when the 
distance between the layers, d, is larger than the inter- 
electron distance so that the discrete nature of the elec- 
trons is irrelevant in the calculation of the interaction 
energy between electrons in different layers |2(j. More- 
over, we restrict ourselves to the case in which the "pas- 
sive" layer has a large electron density, np (i.e. small 
r s ) 1 while the "active" layer has a low electron density, 
7ia (i.e. r s ~ r c , where r c is the critical value of r s for 
the liquid-crystal transition - for a single band with an 
isotropic effective mass, r c is estimated 13] to be r c ss 38.) 
In this limit, the passive layer acts as a ground-plane, so 
the sequence of microemulsion phases in the active layer 
is the same as was derived previously 0,13 for the 2DEG 
in a MOSFET. 

For zero quenched disorder, the FL to WC transition 
is always expected to be first order |14J, In the absence 
of long-range interactions, this would lead to a regime 
of density exhibiting two-phase coexistence; due to the 
screening by electrons in the passive layer, the electrons 
in the active layer interact via a dipolar interaction. It 
has been shown in Q that in this case, first order transi- 
tions are forbidden because the energy of a long interface 
between the two phases is negative. As a result, there 



exists between the WC (for ua < nyyc) an( i the FL (for 
n a > n Fh)i an intermediate range (j^wc < n A < n FL) 
of densities in which there is a set of new "microemul- 
sion" phases which on the mean field level can be charac- 
terized as a mixture of microphase separated regions of 
liquid and crystal [3,13 ■ The relative concentration, sizes, 
shapes, and organization of these regions is determined 
by thermodynamics. Even the full enumeration and clas- 
sification of the distinct phases that result from Coulomb 
frustrated phase separation is not complete, much less a 
thorough investigation of their properties. 

To be concrete, we will focus on the drag effect in 
the case when the active layer is in the bubble phase 
in which the majority phase is a FL, with a finite con- 
centration of bubbles of WC. This phase can be viewed 
as a suspension of pieces of Wigner crystallites floating 
in an otherwise uniform Fermi fluid, like ice in a river. 
(Many of the same considerations apply to more general 
micro-emulsion phases, as we will discuss in a forthcom- 
ing paper [Tj-) As a function of decreasing, i%a, the areal 
fraction of Wigner crystal, fwc rises continuously from 
vanishingly small when ua — n c FL to fwc — 1 when 
nA = n wc . 

The WC bubble phase occurs when ua is close to n FL 
(See Figs.l and 2 in 0] and q)- As ha ^ n FL> the 
bubble size approaches a constant value Lp = dexp[y], 
where 7 is a (positive and order 1) dimensionless ratio of 
microscopic parameters which is proportional to the sur- 
face tension between the Wigner crystal and Fermi fluid 
phases. The spacing between bubbles [n FL — n^J^ 1 / 2 
diverges as ua — > n FL . At mean-field level, these bub- 
bles, themselves, form a crystal of bubbles, but for suffi- 
ciently dilute bubbles, either quantum and classical fluc- 
tuations inevitably melt the bubble crystal. Nevertheless 
the charge inhomogeneities in the resulting bubble liquid 
are still very slowly fluctuating compared to the relax- 
ation times of the Fermi liquid. The competition between 
the local tendency to phase separation and the Coulomb 
interaction (capacitive energy) , results in a mean density 
difference between the WC and FL regions An ~ y/riA/d. 

The qualitative temperature and magnetic field depen- 
dence of fwc is determined by the Pomeranchuk effect 
0, l3 • Because the Wigner crystal has higher spin en- 
tropy density (Swc ~ n l°g[2]) than the Fermi liquid 
(S ~ n(T/Ep)), fwc is an increasing function of tem- 
perature at low temperatures. This is a precise analogue 
of the Pomeranchuk effect in He 3 . Similarly, since the 
spins are substantially more polarizable in the Wigner 
crystal phase, fwc is an increasing function of Hu , but 
it saturates above a (temperature dependent) character- 
istic magnetic field strength, H* at which the spins are 
fully polarized. Moreover, the two effects compete, so 
that when H\\ ^> T, the temperature dependence of fwc 
is quenched. (Here, and henceforth, we will adopt units 
in which h — kp — the Bohr magneton = 1.) 

The strong dependence of fwc on ua, T and Hu is the 
origin of the strong dependences of all other properties of 
the system on these physical variables. To get a feeling 



for the expected dependences, consider the case in which 
T and H\\ are small enough that the Fermi liquid free 
energy can be well approximated by its zero temperature 
value, but large enough that we can ignore the (exponen- 
tially small 15, 16] in powers of -y/rj) magnetic exchange 
energy in the Wigner crystal. In this case, and in the 
absence of disorder, the difference in free energy per unit 
area between the uniform WC and FL phases is 

AF(n A ,r,if||) = AE(n A ) - n A T \og[2cosh(H ]l /2T)}. 

(3) 
Since fwc is determined by the competition between the 
local tendency to phase separation and the Coulomb in- 
teraction, it is a smooth function of AF; we can obtain 
the T dependence of fwc from the T dependence of AF 
by the chain rule. For H» = 0, since AF is a linear 
function of T, it follows that fwc is as well. 

Drag in a dilute bubble liquid: The reasons that 
Rp is so small in a Fermi liquid are not hard to fathom. 
The fluctuations in the charge density in the active layer 
are small when T/Ep is small, and are further reduced 
when averaged over a large length scale, d ^S> 1/y/nX, 
hence the dependances exhibited in Eq. [21 In a mi- 
croemulsion WC bubble phase, by contrast, there are 
large amplitude charge inhomogeneities on the length 
scale, Lb , which implies strong coupling to the electrons 
in the passive layer. 

In the regime when the bubble concentration is small 
there are two types of the current carriers in the active 
layer: electrons and crystalline bubbles. At zero temper- 
ature the electron Fermi-liquid contribution to the drag 
resistance vanishes and the drag effect is entirely due to 
motion of Wigner crystallites. 

Each Wigner crystal bubble in the active layer casts 
an image potential in the passive layer. As bubbles in 
the active layer move with respect to the passive one, the 
electrons in the passive layer scatter on the moving image 
potential. To compute the drag resistance, we compute 
the mean force per unit area, Fpa, exerted on the elec- 
tron fluid in the passive layer when a current J is passed 
through the active layer. In the linear response regime, 
it is clear that Fpa oc ubVb where ns = fwc/^L 2 B is 
the concentration of bubbles and vb is the mean drift 
velocity of the bubbles. Thus 



Ri 



1 A nB 
-~ A D jJLB 

e z riA 



(4) 



where /ib = vb/v is the relative mobility of the bubbles, 
v = J/eriA is the mean electron drift velocity in the ac- 
tive layer, and Ap is a dimensionless constant (discussed 
below) which depends on the properties of the electron 
liquid in the passive layer. 

In the absence of disorder, /.is —> 1- Moreover, as we 
will discuss below, \xb remains non-zero in the presence 
of weak disorder so long as we are dealing with a bubble 
liquid phase. Thus, Rp is similarly non vanishing as 
T — > 0! Moreover, at low temperatures, the dominant T 
and Hu dependances of Rp are through the dependence 



of fwc on these variables - i. e. the Pomeranchuk effect 
discussed above. 

The physics of the proportionality constant Ap is in 
general, complicated - similar complexities occur as in 
the closely related problem of electro-migration in met- 
als. To compute Ap we need to assess the nature of 
the scattering potential induced in the passive layer by 
a bubble in the active layer. The long range tails of po- 
tential are screened away by the electrons in the passive 
layer. The magnitude of the potential, Vb ~ edAn, is 
determined by the magnitude of the charge deficit in the 
active layer, An, associated with the bubble. From the 
previously mentioned scaling of An with d, it follows that 
the WKB barrier penetration factor, W = ^/meVBLB is 
larger than 1, so the potential behaves like a hard object 
of radius ~ Lb being dragged through the passive layer - 
electron tunnelling through the image potential is negli- 
gible. This is one of the key reasons that a bubble phase 
produces a large drag resistance. 

It requires additional microscopic analysis to obtain 
an estimate of Ap, and there are numerous "regimes" 
depending on the relative values of n B (the mean spac- 
ing between bubbles), Lb, and the elastic and inelastic 
mean- free paths, £ e i and £ e - e , in the passive layer. The 
important general point is that Ap is not generally small. 

To give a concrete example, consider the case in which 



,-1/2 



3> l e i S> Lb and £ ee ^> £ e i . In this case, each bub- 



ble acts as an uncorrelated scattering center with cross 
section proportional to Lb, and consequently 



A T 1/2 

Ap = L B n P ' 



(5) 



In a future paper, we will report an analysis of Ap in the 
various other regimes. 

Bubble mobility. In the absence of quenched disor- 
der, there are two distinct bubble phases - the bubble 
crystal and the bubble fluid. 

In the bubble liquid in the absence of pinning, the bub- 
bles are swept along with the fluid, so that fiB = 1 and 
R = Rp. In the presence of weak disorder, we still ex- 
pect that «s sw 1 at temperatures large compared to the 
pinning potential, although in this case, the additional 
channels for momentum relaxation in the active layer will 
generally result in R > Rp . 

The T — > behavior of u B is a more subtle issue. At 
riAd 2 ^> 1 the bubbles are large compared to the spacing 
between electrons. To first order they can be treated 
as classical objects - at this level of approximation, they 
are pinned at favorable sites of the disorder potential. 
Quantum mechanically, however, the bubbles can tunnel 
from one pinning site to another. Moreover, the fact 
that the number of bubbles is only conserved on average 
affects the nature of this tunnelling process profoundly. 
Specifically, since the Fermi liquid degrees of freedom are 
"fast," we can integrate them out to obtain an effective 
action for the bubbles. (This can be done simply, using 
the method of Levitov and Shvtov|l8|. as we will show 
in a forthcoming paper.) Among other interactions that 



are generated by doing this is an effective bubble hopping 
matrix element, t(R,R') ~ \R - R'\~ p where (3 oc 1/G 
and G is the conductance of the Fermi fluid in units of e 2 . 
Physically, this derives from a process in which a bubble 
at pinning site R virtually melts, the associated electron 
density deficit propagates through the Fermi fluid, and 
then recrystallizes at site R'. It is well known [ljj that for 
/3 < D, where D = 2 is the spatial dimension, Anderson 
localization does not occur. This establishes that, for 
clean enough systems, the bubbles, although large, are 
not localized, /is is typically a fairly strongly increasing 
function of T, but it approaches a non-zero value as T — > 

if G is large enough! 

Qualitative comparison with experiment: The 
theoretical model sketched above contains the ingredi- 
ents necessary for a qualitative understanding of much 
concerning the listed experiments: 

The drag resistance in the WC bubble phase Rd is 
not parametrically small, either in powers of 1/kpd or of 
T/Ep, as in the Fermi liquid, which accounts for point 

1 in the above. According to Eq. 5 Rd is propor- 
tional to ub and /xg. Since Lq is T and Hu indepen- 
dent, riB oc fwc, so the general features of the T and 
Hu dependances of Rd and R, points 2, 3, 4, 6, and 
7, follow directly from the Pomeranchuk effect, and the 
interplay between the T and Hu dependence is readily 
understood from Eq. [3J In particular, fwc and, conse- 
quently, Rd are increasing function of Hu which saturate 
at H\\ > H*. When the spins are polarized by H\\ > T, 
they no longer have any entropy, so the T dependence 
of fwc an d the corresponding contribution to the In- 
dependence of Rd (T) are quenched. 

In the case of very small disorder (is is close to unity 
and does not exhibit strong T and H\\ dependances. Then 
R(T,Hu) is entirely determined by the aforementioned 
dependances of tib- The more detailed T-dependances 
of Rd(T, Hu), points 2 , cannot be discussed without a 
microscopic calculation of fiB (T, Hu). . 



Since the dominant T and H\\ dependances of Rd and 
R derive from their implicit dependence of fwc, it is 
unsurprising that these dependances are similar. 

It is an unavoidable consequence of the existence of mi- 
croemulsion phases that, in the absence of quenched dis- 
order, there is a sequence of continuous phase transitions 
(as opposed to a putative first order WC to FL tran- 
sition) between different microemulsion phases and be- 
tween these phases and the low and high density uniform 
phases. These phases become increasingly less conduct- 
ing with decreasing ua, until the limiting WC phase is 
insulating. Thus, for weak enough disorder, an apparent 
metal insulator transition is inevitable, thus accounting 
for the observation that gave the subject its name, point 
6. Since \xb does not vanish as T — + the conceptually 
key conclusion of our analysis is that Rd does not van- 
ish as T -v in the WC bubble phase: R D (T = 0) ^ 0. 
This fact has not been confirmed experimentally, perhaps 
because Rd(T = 0) is small. 

Finally we would like to compare our approach with 
that of Ref. [21| . to the experimental data on the drag 
resistance. The approach of Ref. (2J is supposed to work 
well at high electron concentration far from the regime 
of phase separation. Extending the results of |21| to the 
regime of low n one, in principle, can explain the large 
value of Rd- On the other hand we suspect that the 
dramatic increase of Rd(Hu) in the magnetic field and 
quenching it's temperature dependence in low n samples 
can not be simply explained using that approach. 
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